UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADPO 13247 

TITLE:  Ballistic  Conductance  of  a Quantum  Sphere 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 
Availability:  Hard  copy  only. 


This  paper  is  part  of  the  following  report: 

TITLE:  Nanostructures:  Physics  and  Technology  International  Symposium 
[9th],  St.  Petersburg,  Russia,  June  18-22,  2001  Proceedings 

To  order  the  complete  compilation  report,  use:  ADA408025 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADP013 147  thru  ADPO  13308 


UNCLASSIFIED 


9th  Int.  Symp.  "Nanostructures:  Physics  and  Technology’ 
St  Petersburg,  Russia,  June  18-22, 2001 
© 2001  Ioffe  Institute 


GPLDS.03p 


Ballistic  conductance  of  a quantum  sphere 

J.  Bruningt,  V.  A.  Geyleri,  V A.  Margulis%  and  M.  A.  Pyataevi 
f Institut  fur  Mathematik,  Humboldt-Universitat  zu  Berlin, 

Unter  den  Linden  6,  D- 10099  Berlin,  Germany 
J Institute  of  Physics  and  Chemistry,  Mordovian  State  University, 

430000,  Saransk,  Russia 

Abstract.  The  conductance  of  a quantum  sphere  with  two  one-dimensional  wires  attached  to 
it  is  investigated.  An  explicit  form  for  the  conductance  as  a function  of  the  chemical  potential  is 
found  from  the  first  principles.  The  fonn  and  positions  of  the  resonance  maxima  on  the  plot  of  the 
conductance  are  studied. 


Introduction 

The  aim  of  this  paper  is  the  theoretical  study  of  the  ballistic  electron  transport  of  a nanode- 
vice consisting  of  a sphere  and  two  wires  attached  to  the  sphere.  We  consider  an  idealized 
model  in  which  the  wires  are  taken  to  be  one-dimensional.  This  crucial  simplification  is 
based  on  the  possibility  to  describe  the  electron  motion  in  nanowires  only  by  means  of  lon- 
gitudinal part  of  wave  functions.  To  join  the  wave  functions  in  the  wires  with  those  in  the 
sphere  we  use  boundary  conditions  at  points  of  gluing  the  wires  to  the  sphere.  These  bound- 
ary conditions  are  similar  to  those  in  the  zero-range  potential  theory  and  lead  to  appearance 
of  phenomenological  parameters  like  the  scattering  length  for  a zero-range  potential  [ ].  A 
useful  mathematical  formalization  of  the  approach  considered  here  is  founded  on  the  Krein 
resolvent  formula  from  the  self-adjoint  operator  theory  and  gives  the  scattering  matrix  in 
terms  of  the  renormalized  Green  functions  for  the  free  Hamiltonians  on  the  sphere  and  in 
the  wires. 

1.  Boundary  conditions  and  scattering  matrix 

We  consider  a wire  R^  ( / = 1, 2)  as  the  semi-axis  x > 0.  The  wire  Rj  is  attached  to 
the  sphere  S (S  has  the  radius  R)  by  gluing  the  point  0 from  R+  to  a point  qj  from  S. 
Denote  by  Ho  the  electron  Hamiltonian  on  S,  which  coincides  with  (2m* R2)~lL2,  where 
L is  the  angular  momentum  operator  and  m*  is  the  effective  mass.  A wave  function  / of 
the  device  consists  of  three  parts:  /o,  /i , fi,  where  fo  is  a function  on  S,  and  fj  (j  — 1,2) 
is  a function  on  R~.  Let  r(q,  q')  be  the  geodesic  distance  between  points  q and  q'  on  the 
sphere  S.  Since  we  are  going  to  use  boundary  conditions  of  the  zero-range  potential  theory, 
we  must  consider  the  function  fo  having  the  following  asymptotics  near  the  points  qj  [ ]: 

fo{x)  = ajifo)  \nr(x,qj)  + bj(fo)  + o(l),  (1) 

where  a}  (fo)  and  bj  (fo)  are  complex  numbers;  they  play  the  role  of  the  boundary  values 
for  the  function  fo . As  to  the  functions  fj,  the  role  of  their  boundary  values  are  played,  as 
usual,  by  fj  (0)  and  /j  (0).  For  reasons  of  space  we  shall  consider  here  only  the  boundary 
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conditions  of  the  Ne  umann  type;  the  most  general  form  of  these  conditions  is  given  by  the 
expressions 

bj(fo)=  E [fijkUk(fa)  + Ujkfk(O)] , 

2=1  (2) 
/•( 0)  = E [ajkOkifo)  + Yjkfm]  ■ 

k= 1 

Here  the  parameters  ajk,  fijk  and  yjk  forms  2 x 2-matrices  A,  B,  and  C,  respectively,  such 
that  the  4 x 4 -matrix 


is  Hermitian.  Let  Go(q , q'\  E)  be  the  renormalized  Green  function  of  Ho,  i.e.,  the  function 
obtained  from  the  Green  function  Go(q,  q'\E)  by  extracting  the  singular  term 
(m* R2 / 71  h2)  \nr(q,  q' , ).  Denote  by  Qo(E)  the  Krein  Q-matrix  [ ];  this  is  a matrix  with 
the  elements  QJ0k (E)  = Goiqj , qk\  E),  if  j ^ k,  and  £>jf  (£)  = Goiqj,  qk',  E),  if  j = k. 
Using  the  explicit  form  of  Go  we  can  show  that 

Qo(E)=  Ql2(E)  = -^  Q + t(£))- |tg(trt(£))-ln(2i?)  + C£  . (3) 

where  is  the  logarithmic  derivative  of  the  T -function,  Cg  is  the  Euler  constant,  and 
t(E)  = (2 Vfi2  + 2 m*R2E.  The  non-diagonal  elements  of  Qo(E)  have  the  form 

ei2(^)  = (^)  = (-  . (4) 

where  Vv(x)  is  the  Legendre  function  and  r 12  = r(q\,q2).  The  transition  amplitudes 
Sjk(E)  from  the  channel  Rt  to  K at  the  energy  level  E and  reflection  amplitudes  Sjj  (E) 

in  the  channel  Rt  forms  the  scattering  matrix  on  the  sphere  S(E) . Using  the  Krein  resolvent 
formula  we  can  find  as  S(E): 

~ ifi2  "1  ["  i h2  -1  “ 1 

5=  — k + C + A*  (Q0(E)  - B)~]  A —k-C-A*(Q0(E)-B)-1A  (5) 

_ rrr  J \_m* 

(here  k — \j2m*  E /h  is  the  electron  wave  vector).  Eq.  (5)  shows  that  the  matrix  B is  just 
a matrix  of  parameters  for  a zero-range  perturbation  of  Ho  at  the  points  q/ ; therefore,  to 
avoid  effects  of  non-locality  we  must  choose  B in  a diagonal  form:  fijk  — SjkPj.  It  is 
known  that  the  diagonal  element  Pj  is  expressed  in  terms  of  the  scattering  length  X(-)]  on  the 

zero-range  potential  at  the  point  qj  as  follows:  f$j  — —m*  \n(Ap  )/jif)2.  The  matrix  C has 
a similar  meaning,  namely,  C is  a matrix  of  parameters  for  a zero-range  pert  urbation  of  the 
free  Hamiltonians  in  the  wires.  Thus,  we  must  suppose  C to  be  diagonal  with  the  diagonal 
elements  yj  — — m*X{p /h 2 where  is  the  scattering  length  for  the  zero-range  potential 

at  the  point  0 in  the  wire  . As  to  the  matrix  A,  it  is  responsible  for  the  transmission  from 

the  wires  to  the  sphere.  Indeed, ifA  = 0,  then  the  transition  coefficient  IS12I2  = IS21I2  = 0. 
Eqs.  (2)  show  that  in  the  case  of  ajk  / 0 for  j / k there  are  non-trivial  boundary  conditions 
which  connect  the  wire  RJ  with  the  point  £/*;  therefore,  we  must  suppose  A to  be  diagonal, 
too.  Moreover,  the  matrix  S(E)  is  symmetric  only  in  the  case  of  real  elements  ctjj  = o/j. 
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It  is  convenient  to  express  these  elements  in  terms  of  parameters  Xj  having  the  dimension 
of  length:  a2  — ( m*)2Xj/H 2 . Thus,  in  our  model  the  scattering  on  the  sphere  is  described 
by  means  of  the  six  real  parameters  a j,  fij,  yj.  Note,  that  in  the  case  of  a two-dimensional 
system,  the  zero-range  perturbation  vanishes  in  the  limit  ft  j —*■  oo;  Eq.  (5)  shows  that 
in  this  limit  S\2(E)  — *•  0,  as  might  be  expected.  We  stress  that  the  idealization  of  the 
one-dimensional  wire  is  possible  only  in  the  case  that  the  cross-section  of  the  real  wire  is 
much  less  than  the  typical  sizes  of  the  system.  In  particular,  our  model  works  only  in  the 
case  of  relatively  large  distance  r\i\  namely,  we  shall  suppose  r\ 2 is  vastly  larger  than  the 
Fermi  wave  length  of  an  electron  in  the  wire.  Nevertheless,  we  can  get  the  proper  limit 
(|5i2(£')|  — * ► 1 as  r\2  —*■  0)  using  more  general  boundary  conditions  than  those  in  (2). 


2.  Results  and  discussion 

For  reasons  of  symmetry  we  shall  suppose  y\  = y2  = y and  denote  X ^ simply  by  A.  In 
this  case  Eq.  (5)  gives  the  transition  coefficient  T\2(E)  = | .S’  1 2 ( E ) | 2 in  the  following  form: 


Tn{E)  = 


(\6k)2X\X2  |Qq2|2 

4 X\X2k2  - 2ik(4  - ikX)(X2Q\x  + X\ Qf) 


(4  — ikX)2  det  Q 2 


(6) 


where  Q(E)  is  the  dimensionless  Q-matrix:  Q(E)  = (fi? / m*)(Q(E)  — B).  Using  Eq.  (6) 
we  can  find  the  conductance  G as  a function  of  the  chemical  potential  11  at  temperature  T. 
Namely,  according  to  the  Landauer-Biittiker  formula 


2 00 

G{H,  T)  — J Tn(E)  ) <1E,  (7) 

0 

where  fo  is  the  Fermi  distribution  function. 


Fig.  1.  Transition  coefficient  Tu  as  a function  of  kX  ( k = y/2m*E/h ):  (a)  case  of  a generic  position 
for  wires  (ri2  = ~Jn  R)',  (b)  case  of  an  antipodal  position  for  wires  (ri2  = n R). 
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In  the  discussion,  we  restrict  ourselves  to  the  case  of  Xj  — X(j)}  — X.  The  quantity  T\2 
as  a function  of  E has  a series  of  shatp  splash-like  maxima,  their  positions  are  determined 
by  the  zeroes  of  det  Q(E),  i.e.,  by  the  levels  of  the  zero-range  perturbation  of  Ho  with 
parameters  fij.  On  the  other  hand,  if  r\ 2 ^ : tR,  then  det  Q(E)  has  poles  of  the  second 
order  at  the  points  of  the  spectrum  of  Ho,  i.e.,  at  the  energy  levels  £7  = h2l(l  + l)/2  m*R2 
of  an  electron  on  the  sphere  S,  and  Eq.  (7)  shows  that  T^iE/)  = 0 (see  Figure  1(a)).  It 
follows  from  the  numerical  analysis  that  £12  depends  only  slightly  on  the  parameters  ji 
and  y. 

If  the  points  q\  and  q 2 are  antipodal  (ri2  = it  R),  then  the  behavior  of  £12  at  points  £/ 
changes  drastically.  Indeed,  at  E — E\  the  numerator  and  the  denominator  in  Eq.  (7)  have 
a pole  of  the  same  order,  hence,  T\i  (£7)  does  not  vanish.  Therefore,  the  oscillation  minima 
are  not  positioned  on  abscissa  but  they  lay  on  a curve  of  the  form  7)2  = k(ak2  + bk  + c)-1 
(Figure  1(b)).  In  the  region  / 1 the  oscillation  period  A k is  practically  constant:  A k = 

. A numerical  analysis  of  the  conductivity  at  finite  temperature  shows  that  with  a rise 
of  temperat  ure  the  peaks  are  smoothed  and  their  amplit  udes  are  lowered. 
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